A lattice calculation of the pion form factor with Ginsparg- Wilson-type fermions 
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Results for Monte Carlo calculations of the electromagnetic vector and scalar form factors of the 
pion in a quenched simulation are presented. We work with two different lattice volumes up to a 
spatial size of 2.4 fm at a lattice spacing of 0.148 fm. The pion form factors in the space-like region 
are determined for pion masses down to 340 MeV. 
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I. INTRODUCTION 

The ab-initio determination of hadronic properties is 
an important task of lattice QCD calculations. Follow- 
ing the impressive progress in the determination of the 
hadron mass spectrum, one went on to more advanced 
calculations, such as the determination of matrix ele- 
ments. An important example is the electromagnetic 
pion form factor 0, S i, H SB • 

With the successful implementation of chiral symme- 
try, based on the Ginsparg- Wilson (GW) relation 8] , also 
results at smaller quark masses are feasible. In addition 
to an exact implementation |^ 0| , further Dirac opera- 
tors have been suggested that fulfill the GW condition 
in an ap prox imate way, among them the domain-wall 
fermions [ill , the perfect fermions and the chi- 
rally improved (CI) fermions 0, ^| . 

Dynamical calculations at small quark masses are cer- 
tainly the ultimate goal. Currently such calculations are 
still prohibitively expensive but one can try to improve in 
two ways: either one works with dynamical fermions in a 
conventional formulation or one tries to approach smaller 
quark masses with GW type operators in the quenched 
approximation. In our present contribution we focus on 
the latter. 

We present here a first-principles calculation for off- 
forward lattice matrix elements of operators that mea- 
sure the vector and scalar form factors of the pion. We 
work with the CI fermions, which constitute an approx- 
imate solution of the Ginsparg- Wilson relation and thus 
have improved chiral properties (16| . In our calculation 
we consider comparatively low pion masses down to 340 
MeV, which is likely in the domain where chiral perturba- 
tion theory (ChPT) can be applied. Another convenient 
feature of chirally improved fermions is that the renor- 
malization of the local current Zy is quite close to 1 ^3 ■ 

Preliminary results of our calculations have been pre- 
sented in Ref. 18] . The quenched results for the vector 



form factor at the available transferred momenta lie close 
to the experimental curves. From our data for the form 
factors we can derive the charge radius of the pion as 
well as its scalar radius. The former comes out smaller 
than experiment, due to the too large value of the p- 
meson mass. The scalar radius appears to be sensitive 
to quenching and possibly to the contributions of the 
disconnected diagrams, which we do not include in our 
calculations. 

The pion form factor in the quenched approximation 
has been studied for Wilson type fermions [M S S IE j 
for GW type fermions and for twisted mass fermions 
There are few results for dynamical fermions: for 
Wilson- type fermions in Refs. and a study with 

domain-wall valence quarks and dynamical (MILC) Asq- 
tad fermions P| (and further references therein). 



II. PION FORM FACTORS 



The vector form factor is defined by 
(7r+(p/)|F^K+(p,))co„t = {Pf+P^)^.F^iQ^) , (1) 

where = (pf — pi)^ = —t is the space-like invariant 
momentum transfer squared, and 
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is the vector current. 

The pion vector form factor is an analytic function 
of t. Its space-like values are determined from the val- 
ues on the boundary of its analyticity domain, i.e., the 
cut along the positive real t-axis. Although it starts at 
t = 4mg, significant contributions only come from the 
hadronic region, starting with the two-pion threshold at 
4toJ. Along the first part of the cut, until inelastic chan- 
nels become important, due to unitarity the phase shift 
is essentially the p wave phase shift of the elastic two- 
pion channel. That region is dominated by the p-meson 
resonance. This feature has been exploited by various 
dispersion relation representations _21,, 22, 23] . 

Although there are further inelastic contributions from 
the four-pion channel, these remain tiny until the nuj 
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channel opens. There is also a small contribution from 
the isoscalar w, coupling through higher-order electro- 
magnetic interactions. All these contributions show small 
effects on the near space-like region and can be safely ne- 
glected within the accuracy attained in our work. 

The dominant contribution of the vector meson reso- 
nances has inspired the vector meson dominance (VMD) 
model, where the form factor in the space-like domain is 
approximated by a sum of poles in the time-like region, 



q=p-r 




FIG. 1: Schematic diagram for the matrix element @. 
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In first approximation the coefficients are related to the 
coupling fv-n-TT of the vector meson to the two-pion state 
and its coupling to the photon, fy-f- 

In QCD the large behavior in leading order '23| is 
given by 



8^«,(g2)/2 



(4) 



with as the running coupling constant (logarithmic in the 
argument) and /^r the pion decay constant. In the Q^- 
region accessible to us the logarithm is not identifiable 
and the VMD behavior provides a good approximation. 

Due to electric charge conservation one has F-^ (0) = 1 . 
The mean charge radius squared is defined through 



^ {r\ = 6dF^{t)/dt\t=o 



(5) 



The current PDG average for its value is 0.45(1) fm^ l25|. 

In the simplest VMD model with just the leading p- 
resonance one has (due to normalization) /pTrjr = fpj and 



F4Q^ 



^1/ {ml - t) 
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(r')v,VMD = 0.39 fm" 



(6) 



The scalar form factor is given by the matrix ele- 
ment of the scalar operator, i.e., 

(Tr+{pf)\muUu + mddd\Tr+{p^)) = r^{Q^) . (7) 

Within chiral perturbation theory the scalar form factor 
at = is the so-called sigma term which behaves like 
r7r(0) ^ near zero momentum transfer. The scalar 
radius squared can be obtained from 



(8) 



For a detailed discussion of chiral perturbation theory in 
this context see US- 



UI. STRATEGY 

In order to compute the form factors, we need to eval- 
uate off-forward matrix elements (at several transferred 
momenta), i.e., the expectation values of 



Tr 



^ e'^yS{0,y)OY^ e-'P''5(y,x)75^(x,0)75 
.y;yo=T x;xo=t 

(9) 

where S{y, x) is the quark propagator from x to y and O 
denotes the operator inserted at y. We use the notation 
p = p/, r = Pi and denote the momentum transfer as 
q = p - r (cf. Fig.nj. 

The necessary matrix elements have been evaluated by 
using the sequential source method. The matrix element 
is then written as 

Tr ^ 6^1^5(0,2/) OE(y, 0)75 ■ (10) 
-y;yo=r 

The sequential propagator is defined as 

^-^^P--S{y,x)j,S{x,0), (11) 



E(y,0)= J2 ^ 

x;XQ—t 



and can be easily computed by an additional inversion of 
the Dirac operator D for each choice of the final momen- 
tum p. 



^i?(z,y)S(y,0) = e-^P =' 75 5(^,0) 



Zo=t 



(12) 



Changing the properties of the sink requires the compu- 
tation of new sequential propagators, and so simulating 
several final momenta, different field interpolators, or a 
different smearing for the sink rapidly becomes rather ex- 
pensive. For this reason we have limited ourselves to only 
one value of the final momentum. An alternative version 
of the sequential source method, which builds the sequen- 
tial propagator from the two propagators that sandwich 
the operator, would allow more values of the final mo- 
mentum for free, but on the other hand it would require 
a new inversion of the Dirac operator for every new type 
of operator and for each different value of the momentum 
transfer. This is much less convenient in our situation be- 
cause we consider both the vector and scalar currents and 
want to obtain the corresponding form factors for several 
values of the momentum transfer. 
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We extract the physical matrix elements by computing ratios of 3-point and 2-point correlators: 



Rit, t; p, q) 



(Pft;p)O(r^q)F(0;r)) 
(P(t;p)P(0;p)) 



' {P{t; p) P(0; p)) (F(r; p)_P(0; p)) {Pjt - r; r)P(0; r)) 
{P{t; r) P(0; r)) (P(t; r) P(0; r)) (P(t - r; p) P(0; p)) ' 



r 



(13) 



TABLE I: Values of the non-zero momentum transfers. 



size 


V2po 


2po 




2^2 po 


16^ X 32 
12^ X 24 


0.739 GeV 
0.986 GeV 


1.045 GeV 
1.394 GeV 


1.280 GeV 
1.707 GeV 


1.478 GeV 
1.971 GeV 



where P — uj^d is our pseudoscalar interpolator. We 
keep the sink fixed at time t and vary the timeslice r 
where the operator O sits (scanning a range of times- 
lices) . 

The ratio H13|l eliminates the exponential factors in 
the time variable which are present if one consid- 
ers (P(i; p) 0(t; q) P(0; r)) alone. As a consequence, 
R{t, r; p, q) exhibits two plateaus in r: ^ r ^ t and 
t <^ T <^ T. In the case in which the sink is put at 
t ~ T/2 the ratio is antisymmetric in t — T/2. 

In the Monte Carlo simulations the momentum projec- 
tions of the sink onto non-zero p and the operator onto 
non-zero momentum transfer q lead to significant statis- 
tical fluctuations due to the worse signal-to noise ratio. 
In particular, for many configurations the 2-point corre- 
lators become negative on timeslices near the symmetry 
point t = T/2. A reasonable way out of this inconvenient 
situation is to choose the sink to sit at a timeslice smaller 
than T/2, because in this region the 2-point correlators 
have much larger values, and moreover their relative er- 
rors are smaller. We have thus put the sink at timeslice 
t = 7, while the source remained at timeslice t = Q. The 
matrix elements are finally determined (as in the case 
t = T/2) by combining the two plateau values on either 
side of the sink, Pihs ± Prhs, where the relative sign de- 
pends on the parity properties of the operators. This sign 
is negative for the vector and positive for the scalar form 
factor. 

The term under the square root in Eq. H13|l becomes 
trivial for q = (0, 0, 0). In that case the ratio becomes 



(P(t;p)O(T;q^_0)P(0;r = p)) 
(P(t;p)P(0;p)) 



(14) 



and is sufficient to eliminate the exponentials. 

We choose momenta |p/| = |pi|, which implies Ef — 
Ei so that the transferred 4-momentum is then given by 
= IqP- III this way |3| one achieves for the electro- 
magnetic form factor a cancellation of the kinematical 



factors in ^ and 

(7r(P/)I^Mk(P0)latt 



2^EjE, 



(7r(p/)|y^|7r(p,))c 



(15) 

Indeed, with this choice of momenta, and using the /i = 4 
component of the vector current, the overall factor be- 
comes 



Ef + E^ 

2y/EfEi 



1 



(16) 



This cancellation is unfortunately no longer possible in 
the calculation of the scalar form factor, where a re- 
maining multiplication by the quantity 2 y^Ef Ei is still 
needed. 

In our calculations the pion source is located at the 
origin and the sink has a non-zero 3-momentum which is 
fixed to pf = (l,l,0)po, where po = 27r/(aL) denotes 
the smallest non-vanishing spatial momentum. The mo- 
menta of the current are then chosen such that implicitly 



IP/I 



V2 



Po 



(17) 



is verified. We use twelve different 3-momenta for the 
current, (0,0,0), (0,1, ±1), (1,0, ±1), (2,0,0), (0,2,0), 
(2,1,±1), (1,2,±1), (2,2,0); these are all the possible 
choices that obey the condition IjlTfl . This setup then 
gives rise to four non-zero (and equidistant) values for 
the square of the momentum transfer, = 2npf^ for 
n = 0, 1, 2, 3, 4 (see Table IJ). Momentum conserva- 
tion finally defines the twelve corresponding values for 
the 3-momentum of the source, which all have the same 
module. 

We could also consider larger values of the module of 
the initial momentum, since this would generate further 
values of while still using the same sequential propa- 
gators through the sink. In this case, however, the statis- 
tical errors become much larger, and moreover we could 
not use Eq. (|16|) anymore. A change of the momentum 
at the sink would require the expensive computation of 
a new set of sequential propagators. 



IV. SIMULATION PARAMETERS 

We use the chirally improved Dirac operator 0, 0, 
IT^ . which constitutes an approximate Ginsparg- Wilson 
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TABLE II: Pion and rho masses at the various quark masses 
simulated (taken from 16]). 





16^ X 32 


16^ X 32 


16^ X 32 


12^ X 24 


12=* X 24 


m 


0.02 


0.04 a-^ 


0.06 


0.04 a-^ 


0.06 


m,r 


342 MeV 


471 MeV 


571 MeV 


474 MeV 


575 MeV 


mp 


845 MeV 


895 MeV 


941 MeV 


952 MeV 


976 MeV 



operator. The gauge action is the Liischer-Weisz tadpole 
improved action at /3 = 7.9 [23|, corresponding to a lat- 
tice spacing of 0.148 fm (determined from the Sommer 
parameter in '28'|). Here we do not study the scaling 
behavior. However, although we cannot estimate the 
discretization effects, we expect that our results have 
only small 0(a) contributions, because chirally improved 
fermions are Ginsparg- Wilson- like fermions. This has 
been observed in the study of the hadron spectrum [T^ 
and other hadron properties p9l |. 

The pion interpolators are constructed from Jacobi- 
smeared quark sources and sinks (with k = 0.21 and 
TV = 18), which improves the signal for the ground state. 

We performed simulations on two volumes: 12^ x 24 
(200 configurations for each mass) and 16'^ x 32 (100 con- 
figurations for each mass). This corresponds to spatial 
lattice sizes of a L ~ 1.8 fm and aL ^ 2.4 fm respectively. 

We work at several values of the bare quark mass and 
the corresponding meson masses are given in Table IhI 
(taken from 0|). The jackknife method is used for esti- 
mating the statistical errors of our results. 
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FIG. 2: The ratio 11311 leading to the vector form factor for all 
quark masses and transfer momenta considered (lattice size 
16=* X 32). 



V. RESULTS 

In order to illustrate the quality of the results, we show 
in Figs. |5] and 13 the time dependence of the ratio (|13|l . 
Fitting the central points of the plateaus to a constant 
then leads to the values for the form factors. 

Some combinations of initial and final momenta in- 
duce more statistical noise than others. To increase our 
statistics we average over all combinations that have the 
same momentum transfer. In particular, for = 2pQ 
and — 6pQ we have four times as many data, and for 
Q2 — ipl the averaging doubles our statistics. 

When comparing the results to physical, renormalizcd 



TABLE III: Values of the pion vector form factor (lattice size 



16^ X 32, a = 0.148 fm) obtained for different quark mass 
values. 



Q^[GeV] am=0.02 am=0.04 am=0.06 

1.01(1) 0.996(3) 0.963(13) 

0.546 0.64(15) 0.596(42) 0.580(23) 

1.093 0.54(14) 0.434(42) 0.409(24) 

1.639 0.36(12) 0.296(34) 0.289(20) 

2.186 0.29(14) 0.252(43) 0.239(25) 



quantities we have to multiply with renormalization fac- 
tors relating the raw results to the MS-scheme. For the 
chirally improved Dirac operator these have been deter- 
mined in [13 and all results we show, except for the plots 
of the ratios, are already converted to the MS-scheme. 
We always use values as determined in the chiral limit. 



A. Vector form factor 

Although the vector current is pointlike and not con- 
served, the value Zy = 0.9586(2) turns out to be close to 
unity. This number converts the lattice bare results at 
a = 0.148 fm to the corresponding renormalized contin- 
uum quantities in the MS scheme at the scale of 2 GcV. 
The fact that the renormalization constant Zy is so close 
to 1 means that we have better control over this kind of 
systematic effects. 

The resulting i^7r(0) is not constrained to unit value, 
but comes very close. This fact is obvious from the re- 
sults summarized in Table IIIII These values have been 
computed by combining the plateau mean values as dis- 
cussed before; these have been obtained by averaging the 
values a.t t — 3 and 4 of the l.h.s. plateau and the values 
at t = 20 — 23 of the r.h.s plateau. The errors have been 
determined with the jackknife method. The larger error 
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FIG. 3: Our results for the pion vector form factor F-^ 
compared with results from experiments (lattice size 16^x32). 



FIG. 5: The inverse vector form factor for the three available 
masses on the 16^ X 32 lattice. 
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FIG. 4: Our results for the pion vector form factor com- 
pared with results from experiments (lattice sizes 16^ X 32 and 
12=* X 24). 



bars in some of the plots are an indication of the rather 
poor signal-to-noise ratio for some timeslices close to the 
lattice center. These timeslices are however outside the 
regions that we use for the fits. As already mentioned in 
Sect, mil we profit from the cancellation of kinematical 
factors by following the approach suggested in Ref. 4] . 

Fig. |31 compares our results for the vector form factor 
with results from experiments HSmiHlll. In Fig. H 
we show these results together with those for the small 
lattice size 12^ x 24. For this smaller lattice size the 
statistics was sufficient only for the data at am = 0.06. 
The resulting form factor values are at other momentum 
transfer values than those of the larger volume, but in 
reasonable agreement. 

In Fig. [51 we plot the inverse of the vector form fac- 



tor. The electromagnetic pion form factor in the time- 
like region is dominated by the p-meson. In the near 
space-like region it may therefore be well approximated 
by a monopole form © such that the leading behavior 
of is linear. This is indeed observed in Fig.|Sl 

As the VMD model is just an approximation, one 
expects corrections due to other resonances and more- 
particle channels in the time-like region. In the fit these 
may be taken into account by further pole terms. Our 
data shown in Fig.|Sldo not really require such a multi- 
parameter fit and we therefore discuss only the results of 
a linear fit. 

The derivative of the form factor at = gives the 
charge radius, as shown in Fig.El The derivative has been 
obtained from the linear fit to the inverse form factor. 
Whereas the large mass results are compatible with a 
constant in a m, the number for the lowest mass is smaller 
but has a very large statistical error. These values are 
quite compatible with numbers from Q at comparable 
quark masses. 

Obviously in our results there is an unknown system- 
atic error due to the quenched approximation. We re- 
frain from applying QChPT extrapolations because the 
data are not accurate enough to reliably determine the 
unknown expansion parameters for the non-linear terms. 
Computing the average gives (r^)v = 0.291(18) fm^. 

The value of (r^)v in the VMD model is inversely pro- 
portional to the mass squared of the p-meson. We obtain 
"^VMD — 0.80(5) GeV^. This agrees with the range of 
values obtained for in the analysis of Ref. [I^ in the 
direct p-channel (e.g., to^ = 0.80 GeV^ at am = 0.04). 
These large values for the p-meson mass might at least in 
part explain (via the VMD model) our low (as compared 
to experiment) result for (r^)v. 

Like most other calculations we could afford only to 
work for a single value of the lattice spacing and therefore 
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FIG. 6: The mean square radius for the vector form factor 
as obtained for the three available masses on the 16^ x 32 
lattice. The asterisk denotes the average experimental value 
[25|. the error band shows the result of a constant fit to the 
three points. 



TABLE IV: Values of the pion scalar form factor (lattice size 
16^ X 32, a = 0.148 fm, without disconnected contributions) 
obtained for different quark mass values. 



0'[GeV] 


am=0.02 


am=0.04 


a m=0.06 





1.77(2.20) 


1.30(39) 


1.42(17) 


0.546 


2.09(53) 


1.38(13) 


1.31(6) 


1.093 


1.60(51) 


1.31(13) 


1.24(7) 


1.639 


1.61(56) 


1.13(13) 


1.04(7) 


2.186 


1.37(69) 


1.03(16) 


0.96(9) 
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FIG. 7: The ratio (1131 for the matrix element ||7J leading 
to the scalar form factor for all quark masses and transfer 
momenta considered (lattice size 16^ X 32). 



we have no control on scaling violations. However, due 
to the (approximate) GW nature of the CI operator we 
expect 0{a) corrections to be very small (16|. 

Due to the chosen method we have fewer momentum 
transfer values than who also work at larger lattices 
and at slightly smaller quark masses. Our results, ob- 
tained with different action and Dirac operator, are in 
generally good agreement with the findings of other lat- 
tice calculations in the discussed range of values for mass 
and momentum transfer 0, 0] , even for those with dy- 
namical background Higher quark masses overesti- 
mate the values of the form factor, as it is expected due 
to the corresponding higher vector meson mass. This 
then leads to an underestimation of the charge radius. 



B. Scalar form factor 

The scalar form factor also has disconnected contribu- 
tions which we disregard here (like it is often done due 
to the inherent technical complications of backtracking 
loops). We can notice from the upper plots in Fig. [7| and 
the upper line in Table HVI that the scalar form factor at 
zero momentum transfer divided by the quark mass has 
nearly the same value for all masses that we have simu- 
lated. This means that r7r(0) ^ ~ m, as expected. 



We may determine the scalar radius squared from a fit 
to the Q^-dependence of the scalar form factor (now ex- 
plicitly normalized according to Eq. ©). 

In analogy to the vector form factor we expect that the 
scalar form factor in the space-like region is a decreasing 
but positive definite function. Thus a linear fit to its 
inverse seems appropriate. However, the result for the 
scalar radius squared is quite sensitive to this assump- 
tion. In Fig.|Slwe compare the results of such a fit with 
those from a direct linear fit to our data in the space- 
like region and find a difference of almost 50%. In either 
case, the resulting value may then be extrapolated to the 
chiral limit. 

ChPT relates the pion decay constant to the scalar 
form factor radius via 



with 



4 7r/^ 



(18) 



(19) 



The chiral expansion of the pion decay constant should 
behave like ^] 



fjf=i+^h+o{e) , 



(20) 
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FIG. 8: The mean square radius for the three available 
masses on the 16^ x 32 lattice for the scalar form factor. Up- 
per plot: Determination from the linear fit to the values of 
r,r(<3^)/r,r (0). Lower plot: Determination from the linear fit 
to the inverse of that ratio. 



where the value £4 = — \n{M^ 
trinsic QCD scale A and in 
A « 47r/^; in Ref. -2^ a value of 
Eq. H18|) may be translated to 



r^A^) 
Hit 



) depends on the in- 
is suggested to use 
i 4.0 ±0.6 is quoted. 



(21) 



Recent values quoted in that context are {r'^)s — 
0.61(4) fm^ m or 0.75(7) fm^ 

In QCD one expects correction terms with a logarith- 
mic singularity in the valence quark mass m. As pointed 
out in (36i |. the leading order logarithmic term m logm 
of ChPT involves quark loops that are however absent 
in the quenched case. There will be non-leading loga- 
rithmic terms, though. One therefore could allow a term 
logm in the extrapolating fit. In the range of mass 
values studied here, the statistical accuracy of the data is 
too poor to render such an extrapolation significant. We 
therefore exhibit only the results of a constant extrapo- 
lation. 

The resulting values (Fig. |H|) are almost an order of 



magnitude smaller than for the electromagnetic case. 
From the data in Table Hvl which takes into account the 
renormalization of the scalar operator, Zs ~ 1.1309(9), 
as well as the kinematical factors, we obtain an aver- 



age for the scalar radius of 



0.054(16) fm^ (based 



on the linear fit to the inverse form factor). This is also 
smaller than the values expected for full QCD. Results for 
the pion decay constant in a recent full QCD lattice calcu- 
lation iSli via (CHll lead to the value (r^), = 0.5±0.1 fm^ 
A corresponding analysis of the quenched BGR data |29l | 
gives again a small value, 0.08 — 0.13 fm^. This seems 
to imply that (r^)s is very sensitive to quenching and 
possibly to the omission of the disconnected pieces. 

A recent unquenched computation with Nf — 2 clover 
fermions 19], which uses matrix elements of the scalar 
current and not the pion decay constant, quotes a larger 
value for the scalar radius squared, {r^)s = 0.60(15) fm^. 
However, this number results from a ChPT-motivated 
extrapolation to data that is actually of smaller size and 
for pion masses > 550 MeV. 



VI. CONCLUSIONS 

Chirally improved fermions provide a framework for a 
first-principles lattice study of the hadron structure at 
comparatively small pion masses. In that context we 
present here our investigation of the vector and scalar 
form factors of the pion. We combine the method of 
Q (which for the vector form factor avoids calculation of 
the pion energy and thus removes one source of statistical 
error) with a GW type action, that allows one to reach 
small pion masses for comparatively small lattice size. 

The vector form factor describes the electromagnetic 
structure and our results are generally consistent with 
expectations of its general form. In view of the quenched 
approximation we expect deviations from the experimen- 
tal numbers. The vector charge radius indeed comes out 
40% smaller than in experiment and the scalar radius 
is much lower than one would expect from unquenched 
calculations. 

Further progress could be achieved - as usual - by using 
better statistics, more momenta, and larger lattices, both 
in lattice units (to study finite volume effects) and in 
physical units (to access lower transferred momenta). 
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